Twins and their boundaries during homoepitaxy on Ir(lll) 
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The growth and annealing behavior of strongly twinned homoepit axial films on Ir(lll) has been 
investigated by scanning tunneling microscopy, low energy electron diffraction and surface X-ray 
diffraction. In situ surface X-ray diffraction during and after film growth turned out to be an 
efficient tool for the determination of twin fractions in multilayer films and to uncover the nature 
of side twin boundaries. The annealing of the twin structures is shown to take place in a two step 
process, reducing first the length of the boundaries between differently stacked areas and only then 
the twins themselves. A model for the structure of the side twin boundaries is proposed which is 
consistent with both the scanning tunneling microscopy and surface X-ray diffraction data. 

PACS numbers: 61.72.Nn, 61.72.Mm, 61.72.Cc, 68.37.Ef , 61.05.cp, 61.05.-a 



I. INTRODUCTION 

Stacking faults and twin boundaries are known to se- 
riously affect the properties of thin films and devices. 
The boundaries with the initial crystallite along the twin 
plane is coherent, while the others, namely the side 
boundaries, are usually not. Incoherent boundaries dis- 
turb the thin film rather strongly and may give rise to 
pronounced property changes, e.g. increase of the electri- 
cal resistivity. The magnetic properties of thin films are 
known to crucially rely on the stacking sequence 1 ' 2 . Co- 
herent boundaries may act as electron scattering planes 
in metal nanowires 3 -, as traps for electrons in SiC^£ or to 
limit the performance of optoelectronic III-V nanowire- 
based devices 6 . 

Due to the low energy of their formation, stacking 
faults are among the most frequent defects in thin filmsi 
For the same reason, there is only a small driving force 
for their removal and twins are thus among the most sta- 
ble defects. There is little knowledge on the kinetics of 
stacking fault healing at the atomic level, on tempera- 
ture induced transformations of side boundaries, or on 
twin coarsening. 

Here we use homoepitaxy on Ir(lll) to gain insight 
into the kinetics and atomistics of twin crystallite anneal- 
ing in thin films. In a homoepit axial system no stack- 
ing faults should be present due to energetics. Com- 
plications due to chemical inhomogeneities or epitaxial 
strain are absent. This article is a follow up of previ- 
ous work on the formation and proliferation of stacking 
faults and twin crystallites in thin homoepitaxial Ir films 
on Ir(lll)£— . The present study combines the surface 
view of a partially twinned film obtained by low energy 
electron diffraction (LEED) and scanning tunneling mi- 
croscopy (STM) with the surface+subsurface, quantita- 
tive view provided by surface X-ray diffraction (SXRD). 

We find that the twinning in Ir thin films on Ir(lll) 
gives rise to fractional step heights. Fractional step 



heights have already been observed several times on (111) 
surfaces of face-centered cubic (fee) metals by scanning 
tunneling microscopy^r— . They result e.g. from stack- 
ing fault tetrahedrons, Lomer-Cottrell locks or slip on 
{111} planes inclined with respect to the (111) surface. 
However, for all these cases the fee stacking sequence is 
preserved on both sides of the fractional step, while in 
the present case of homoepitaxial films on Ir(lll) frac- 
tional steps separate twinned and untwinned parts of the 
crystal. 

Twin and stacking fault healing takes place in a two 
step process, in which first the amount of side boundaries 
is reduced through a crystallite coarsening process before 
eventually the amount of twin crystallites decreases. The 
geometry of the side boundaries is characterized by the 
SXRD measurements. 



II. EXPERIMENTAL 

The STM and LEED experiments have been per- 
formed in an UHV system with a base pressure below 
5 x 10 -11 mbar. The sample was cleaned by repeated cy- 
cles of sputtering with 1.5 keV Xe + or Ar + ions at 1100 K 
and annealing to 1600 K. Iridium was deposited from a 
current heated Ir wire with a standard deposition rate of 
1.3 x 10 _2 ML/s, where 1 ML (monolayer) is the surface 
atomic density of Ir(lll). Imaging was performed with a 
home built, magnetically stabilized STMi£. Unless oth- 
erwise specified, all images are differentiated and appear 
as illuminated from the left. The topographs were digi- 
tally post-processed with the WSxM software^. For the 
LEED I/V analysis, a rear view camera LEED system 
was used. 

The SXRD experiments have been performed at the 
SUV instruments^ installed at the BM32 beamline at 
the European Synchrotron Radiation Facility (ESRF), in 
a UHV system, with a base pressure below 3 x 10 -10 mbar 
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coupled to a Z-axis diffractometer. Sample preparation 
and deposition in this UHV system followed the same 
procedure as for the STM and LEED experiments. A 
monochromatic 18keV photon beam incident under an 
angle of 0.273° with respect to the surface was used. The 
corresponding X-ray attenuation length for Ir is 51 A 20 
(23 ML). Note that this is an average value because the 
sample surface had a bowing of 0.2° resulting in a spread 
of incident angle. The incident beam was doubly focussed 
to a size of 0.3 x 0.2 mm 2 (full width at half maximum 
in horizontal and vertical directions respectively) at the 
sample location. Detector slits, located 570 mm away 
from the sample, were set at 2 mm parallel to the sample 
surface (which was vertical), and 2.2mm perpendicular 
to it (with 5 mm guard slits at 200 mm), resulting in an 
angular acceptance of 0.2° in both directions. 

For the general description of planes and directions, 
the standard cubic system is used (a/ r = 3.8392 A). The 
surface plane is indicated by (111) and the corresponding 
direction by [111]; side- facets of the same symmetry but 
not parallel to the surface are indicated by {111} and the 
corresponding directions as (111). 



III. STM EXPERIMENTS 

Fig. DJa)-(d) displays characteristic elements of a ho- 
moepitaxial growth sequence on Ir(lll) already discussed 
in more detail previousl y 1Q i 11 . After deposition of 0.2 ML 
as in Fig. QJa) dendritic single atom high islands are vis- 
ible. They display a triangular envelope due to the pre- 
ferred formation of {lll}-microfaceted step edges. The 
triangular shape of the island in the image center is mir- 
rored (or rotated by 180°) in orientation. This mirrored 
orientation marks stacking fault islands with the atoms 
residing in the threefold coordinated hollow sites corre- 
sponding to an hep rather than an fee stacking sequence. 
Continued deposition gives rise to island coalescence and 
most stacking fault islands switch to regular stacking 
during this process 9 . At the boundary between regu- 
lar and remaining faulted areas narrow stripes of sub- 
atomic width remain, which offer fourfold coordinated 
adsorption sites. This gap and thus also the faulted 
area bounding on one of its sides become stabilized by 
adatoms forming a monatomic width decoration row ori- 
ented along (110) [Fig. 0Jb)]. Continued growth takes 
place preferentially by heterogeneous nucleation at these 
decoration rows. It causes the formation of additional 
faults (fault proliferation^), stabilizes existing ones and 
leads to the embedding of twins. In Fig.[TJc) the situation 
after 10 ML deposited is shown, where overgrown deco- 
ration rows mark the separation of differently stacked 
surface areas. Eventually, further deposition gives rise to 
an irregular rough surface dominated by mounds. They 
result from heterogeneous nucleation at the boundaries 
between differently stacked areas of the films [compare 
Fig. Ed)]. 

The same film is shown in Fig. [TJe) in a demagni- 



fled view. Each of the bright features marks a (110) 
oriented boundary between areas of different stacking 21 . 
In an annealing sequence with succesive annealing in- 
tervals of 180 s this rough and heavily twinned film is 
subsequently heated to higher and higher temperatures. 
In Fig. [TJf ) after annealing to 850 K the film has be- 
come much smoother and the typical structure size has 
increased. Decoration rows or bright features oriented 
along (110) are now absent. After annealing to 1025 K 
[Fig. DJg)] the film is very flat and the characteristic 
structure size is further increased. Straight steps pre- 
cisely oriented along the (110) directions are visible. 
These straight steps have heights of only a fraction of 
a regular monatomic step. In addition, curved steps of 
monatomic height are visible. They are pinned where 
they touch fractional steps. The monatomic steps are in- 
variably curved outward towards the downhill side. The 
outward curvature appears to imply an enhanced surface 
chemical potential, which for the case of Fig. QJg) may 
be traced back to decaying adatom islands also visible. 
The last two annealing steps to 1200 K [Fig. Dili)] and 
1375 K [Fig. I^i)] increase the typical structure size fur- 
ther, i.e. the step density decreases. The distinction 
between rounded monatomic steps and straight steps of 
fractional height is obvious. As indicated by arrows in 
Fig.QJh) and (i) certain areas bounded by fractional steps 
are misoriented with respect to the global [111] orienta- 
tion. This is apparent through twisted steps being partly 
up and partly down steps which indicates the presence of 
dislocations in the film. The fact that even in Fig. DJi), 
after the decay of all adatom islands, pinned monatomic 
steps are invariably curved outwards is possibly caused 
by less strongly bound atoms in the defect structures as- 
sociated with the fractional steps. The defect structures 
define an enhanced chemical potential visible through the 
step curvature. 

Fig. [2fa) shows a height image (not differenti- 
ated) of an area with straight (fractional) and curved 
(monatomic) steps in detail. Example line scans around 
the junction points marked with the arrows in (a) are 
shown in Fig. EJb) and (c). The situation at a typical 
triple junction is depicted in a three dimensional view of 
Fig. Efd). Neglecting for the moment heavily distorted 
areas, only three different kind of steps are found: (i) 
steps with a height of 2.22 A ±0.17 A, i.e steps with a 
height corresponding to the (111) layer distance of 2.22 A; 
(ii) steps with a height of 0.74 A ± 0.09 A i.e. steps 
which display within the limits of error a fractional height 
of 1/3 of a monatomic step; (hi) steps with a height of 
1.49 A ± 0.24 A i.e. steps which display within the lim- 
its of error a fractional height of 2/3 of a monatomic step 
[compare Fig. EJb) and (c)]. The measured step heights 
do not depend on tunneling voltage and are thus a to- 
pographic effect. Evidently, fractional steps cannot be 
present on the (111) surface of a perfect Ir crystal and 
their occurrence proves the presence of extended defect 
structures in the Ir film even after annealing to 1375 K 
[compare Fig.QJi)]. In the vicinity of the fractional steps, 
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FIG. 1: (Color online) (a)-(d) STM topographs illustrating characteristic features at different stages of homoepitaxy on Ir(lll) 
after deposition of (a) 0.2 ML, (b) 0.9 ML, (c) 10 ML, (d) 90 ML at 350 K. Image size in (a)-(d) is 64 nm x 64 nm. (e) larger 
scale STM topograph of a 90 ML film grown at 350 K. (f)-(i) STM topographs of the film shown in (e) after succesive 180 s 
annealing steps at (f) 850 K, (g) 1025 K, (h) 1200 K, (i) 1375 K. Image size in (e)-(i) is 480nmx 480nm. Inset size in (f) is 
120 nm x 120 nm. Small white arrows in (h) and (i) indicate emerging screw dislocation lines (junctions between up and down 
steps) . 




10 20 30 

X(nm) 



FIG. 2: (Color online) (a) Greyscale STM topograph after deposition of 90 ML film at 350 K and additional annealing to 
1025 K for 180 s. Image size is 128 nm x 128 nm. (b), (c) Line scans along paths indicated in (a) illustrating full and fractional 
step heights as well as strain fields associated with fractional steps, (d) Three dimensional visualization of the junction marked 
in (a). 
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FIG. 3: (Color online) STM topographs after deposition of 
90 ML at 350 K, subsequent annealing to 1200 K for 180 s and 
final 0.2 ML Ir deposition at room temperature. The pre- 
dominant orientation of the triangular envelope of the small 
dendritic islands allows one to identify the stacking sequence 
of the underlying film, fee areas are shaded green, faulted ar- 
eas are shaded red. Image sizes are (a) 255 nm x 255 nm and 
(b) 44 nm x 44 nm. 



typically an upward (downward) bending of the upper 
(lower) terrace perpendicular to the step direction takes 
places. It extends over a distance of about 9nm with 
a significant scatter and indicates the presence of strain 
next to the fractional steps. 

A post decoration technique was employed to assess 
the stacking of the surface areas of the Ir films after 
annealing 22 . 0.1-0.2 ML Ir were deposited on the an- 
nealed film surface at room temperature. Dendritic Ir 
islands with the characteristic triangular envelope grow, 
which display the same orientation in an area of uniform 
stacking with a preference of about 90 %. Areas of mir- 
rored stacking display consequently mirrored preferences 
of the triangular island envelopes. Application of the 



post decoration method to the surface after annealing 
at 1200 K is exemplified in Fig. [3l The areas have been 
marked depending on the stacking sequence identified. 
We find that even after annealing to 1200 K more than 
50 % of the surface area of the 90 ML films is twinned 
with respect to the bulk crystal. This result is backed up 
by quantitative LEED I/V analysis 11 and further sup- 
ported in this article by quantitative analysis of SXRD 
measurements which are sensitive to the stacking of the 
film grains over their full thickness. 

As illustrated in Fig.[3jb), fractional steps separate ar- 
eas of different stacking sequence. In fact, for an annealed 
film adjacent areas of different stacking are always sep- 
arated by a fractional step. Thus fractional steps mark 
the lateral border of the differently stacked grains nor- 
mal to the twin plane, i.e. they indicate the positions of 
side twin boundaries. The orientation of the separating 
steps is along the (110) which is also the intersection be- 
tween {111} planes and the (111) surface. As these frac- 
tional steps separate differently stacked grains, the pla- 
nar defects associated with them must be different from 
the cases of fractional step heights analyzed previously 
by STMi^r— . In order to understand the exact nature 
of these boundaries, sub-surface information is required, 
which SXRD provides us with. 



IV. SXRD EXPERIMENTS 

Figure 2] displays the reciprocal lattice of a perfect 
(bulk) fee crystal (green disks), cut in the plane per- 
pendicular to (111) and containing the [211] direction. 
For a truncated crystal, the loss of symmetry in direct 
space creates a reciprocal lattice of crystal truncation 
rods [CTRs, solid lines in Fig. [4] perpendicular to the 
surface and connecting the reciprocal lattice points cor- 
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FIG. 4: (Color online) Reciprocal space cut in the plane de- 
fined by the [111] and [211] directions of an fee crystal. The 
medium grey (green) disks depict the reciprocal lattice of a 
perfect fee crystal; the dark grey (red) disks depict the recip- 
rocal lattice of a twin crystal, i. e. of an fee crystal rotated by 
180° around a [111] axis; the light grey (yellow) disks depict 
the reciprocal lattice of this twinned crystal, but addition- 
ally rotated by 180° around the [511] direction of the initial 
crystal (twin-of-twin) which matches the [111] direction of the 
twin crystal. For multi-twin peaks, half- or quarter-disks are 
employed. Shaded areas represent the truncated crystals: the 
grey crystal is terminated by a (111) surface while the yellow 
one corresponds to the (111) termination of the twin, i.e. a 
(511) surface in the frame of the initial crystal. These trun- 
cated crystals produce CTRs which are displayed as solid lines 
(grey or yellow depending on the crystal). 



responding to the bulk crystal^. Accordingly the main 
intensity maxima coincide with the fee bulk reflections. 
We now consider Ir ho moepit axial growth. As we have 
shown previously^ - — , the overgrown Ir may either pro- 
long the fee stacking, or grow with a stacking fault initi- 
ating the growth of a twinned crystal portion, i.e. a crys- 
tal portion rotated by 180° around a [111] axis. The first 
stacking contributes to the CTRs intensities in a similar 
way as the substrate does, i.e. with intensity maxima for 
the fee bulk reflections. The second stacking gives a re- 
ciprocal lattice rotated by 180° around [111]. As an illus- 
tration, for the CTR passing through 1/3(422) [labeled 
(I) in Fig. g], intensity maxima at (111), (022), (133), 
etc, corresponding to the fee stacking (green disks), are 
supplemented by maxima shifted with 1/3 [111], corre- 
sponding to the twinned fee stacking (red disks). Moni- 
toring the relative fee and twinned intensities along the 
CTRs during growth qualitatively informs on the relative 
fec/twinned fractional the signal of the film parts with 
fee stacking being initially difficult to separate from the 
much higher superimposed signal due to the substrate. 
As extensively shown in the literature for the case of a 
single or a few atomic layers, more careful analysis, no- 
ticeably accounting for interferences between the twinned 
and untwinned regions, allows one to quantitatively de- 



termine the stacking composition in these layers. Below 
we show that such quantitative analysis can also be per- 
formed for much thicker films (^70 atomic layers). 

Figure [5fa-f) show the X-ray scattered intensity along 
a set of six CTRs, all of them along the [111] direction 
[Figs. E(a,c) are the CTRs marked (I) and (II) in Fig. gj 
the others could be sketched in other appropriate cuts 
in reciprocal space], after deposition of 68 ML of Ir at 
350 K. The occurrence and location of strong maxima 
different from the bulk fee reflections along four of the 
CTRs [Figs. [5ja-d)], together with the absence of notice- 
able extra maxima along the two others [Figs.[5fe,f)], are 
an unambiguous signature for a large twin fraction in the 
Ir thin film. In between the Bragg peaks, the intensity 
variations correspond to the roughness only for Fig. [5fe) 
and to a mixture of the roughness contribution and of 
the twin signal for Figs. [5fa,c). The Bragg peaks for the 
substrate [on the grey curve in Figs. [5fa)] are narrower 
than those from the fee regions in the thin films, con- 
firming that the X-ray beam probes prominently the top 
region of the sample where the effective height of the fee 
domains is reduced - as compared to the case for the 
bare substrate - due to the presence of twin domains. 
The appearance of the twin signal was monitored in situ 
during the growth by measuring the integrated intensity 
around the left position marked "twin" in Fig.[5jb). The 
result is shown in Fig. [5fg), where a smooth increase is 
observed starting from zero. Note that in contrast to the 
STM analysis, this evolution is not only taking into ac- 
count the fraction of twins at the topmost layer, but also 
the twin in the layers below as a result of the penetra- 
tion depth of X-rays. A quantitative analysis, presented 
in the next paragraph, will allow us to determine the Ir 
film stacking composition profile and roughness. In ad- 
dition to the typical twin signals, fainter extra features 
are also observed, to more or less extent, along the CTRs. 
These features are marked by vertical arrows in Figs.[5fa- 
d). Some of them are very faint, two of them are even 
barely seen. Their observation, though very sensitive to 
the alignment of the goniometer axis, is confirmed in a 
second series of experiments. Their origin is discussed 
latter. 

The quantitative analysis of the Ir film composition 
was performed by using a simple model which is detailed 
in appendix A, allowing us to simulate the scattered in- 
tensity along the whole set of CTRs with a single set of a 
few parameters that are refined by a least squares mini- 
mization algorithm. Through this procedure we obtained 
the thin film stacking composition in terms of twin/fee 
fractions as a function of the height in the film as dis- 
played in a graphical way in Fig. [5fh). The model finds 
a twinned fraction of 70 % at the surface, which agrees, 
within the limits of uncertainties, with the STM data 
presented latter. We also obtained an rms roughness of 
4.3 A (assumed identical for twinned and fee regions), for 
which a comparison to the data obtained through the 
STM analysis (2.5 A) is hazardous (see Appendix A). 

The large amount of time needed for the measurement 
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FIG. 5: (Color online) (a-f) Integrated intensity scattered along six CTRs perpendicular to the (111) surface, i.e. along the 
[111] direction, after deposition of 68 ML Ir on Ir(lll) at 350 K. Dots represent the data, while solid lines correspond to the 
fit of the whole set of CTRs. (a) and (c) correspond to the CTRs marked (I) and (II) in Fig. [H The Bragg reflections for the 
fee untwinned crystal are indicated in green. Positions for the twinned crystal are indicated by "twin". The vertical arrows 
point to positions where faint extra contributions are observed. The data set displayed in grey in (a) is the integrated scattered 
intensity for the substrate prior to the growth of the Ir thin film, (g) Integrated intensity at the left position marked "twin" 
in (b), as a function of the deposit amount. The dots represent the data, the solid line was computed through the parameters 
from the fit of the CTRs. (h) Fraction of twin/fee regions in the Ir overgrown film, as a function of the layer index, as deduced 
from the fit of the CTRs. 



of a full set of CTRs prevented us from doing such a 
measurement during the growth of the thin film. As an 
alternative we measured the integrated intensity over a 
restricted range in reciprocal space, below a twin peak 
[Fig. [5(g)] and an fee peak as a function of the film thick- 
ness, which a priori only provides us with a qualitative 
characterization of the evolution of the amount of twins. 
This evolution may however be compared to the scattered 
intensity calculated with a model Ir thin film of increasing 
thickness which stacking composition follows the quan- 
titative picture developed as described in the previous 
paragraph. The calculation nicely reproduces the data 
as shown in Fig. [5(g), which further confirms the validity 
of the quantitative picture of the thin film. The evolution 
displayed in Fig. [5(g) may then be interpreted by inspec- 
tion of our model [Fig. [5(h)]. After deposition of about 
10 ML the twin fraction becomes noticeable and increases 
faster and faster. We know from STM experiments^, 
that after 10 ML the growth mode changes from layer- 
by-layer growth to rough growth dominated by heteroge- 



neous nucleation at decoration rows [compare Fig. H(c)] 
with rapid proliferation of stacking faults. 

We then explored the possible stacking sequences on 
the side facets of fee and twin Ir islands, i.e. on their 
{111} facets other than (111) (for twin Ir islands these 
are {115} planes if expressed in the fee frame), and the 
possible appearence of new Ir orientations correspond- 
ing to twinning with these inclined facets as twin planes. 
This leads to six new possible twin-daughter orientations 
of Ir, three from the fee parent and three from the (111)- 
twin parent. Accordingly we explored lines in reciprocal 
space being perpendicular to the {111} and {115} planes 
respectively. Figure [6(a) shows the integrated scattered 
intensity along the line marked (III) in Fig. 2J i.e. along 
the [511] direction which is normal to a {111} plane of a 
twin crystallite. In addition to the expected twin peaks, 
new ones develop as highlighted by the vertical arrows 
and better visible on the grey curve which was obtained 
after annealing. The intensity of such peaks was mon- 
itored during growth [Fig. [6(c)]. These new peaks are 
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FIG. 6: (Color online) (a,b) Intensity scattered along lines in the [511] (a) and [111] (b) directions, after deposition of 68 ML 
Ir on Ir(lll) at 350 K. (a) corresponds to the CTRs marked (III) in Fig. [4] The Bragg reflections for the fee untwinned crystal 
are indicated in light grey (green), those for the twinned crystal are marked by "twin". The vertical arrows point to positions 
of the twins-of- twins (a) and of the side- twins (b), which are better visible for the scans (grey lines, shifted upwards for clarity) 
corresponding the scattered intensity after annealing at 750 K. The intensity measured off the peaks is a mixture between the 
side facets [(511) and (111)] signals and that from the corresponding stacking faults. (c,d) Intensity at the central position 
marked by a vertical arrow in (a),(b) as a function of the Ir deposit (see text). 



generated by twin daugthers from the (lll)-twin parent 
(twin-of-twin) , as can be seen in Fig. HI A similar ob- 
servation is made along a line lying perpendicular to a 
(111) plane [Figs. [6{b,d)]: the fee reflections are supple- 
mented by extra features developing during growth which 
are generated by one of the twin-daughter from the fee 
parent (side- twin), more precisely the one resulting from 
a 180° rotation around the [111] direction. 



V. TEMPERATURE DEPENDENCE 

Figure [7] compares the annealing temperature depen- 
dence of twin fraction estimated by the different tech- 
niques. The black squares in Fig.^a) display the fraction 
of twinned surface area as a function of annealing temper- 
ature measured by LEED. It is around 0.6 up to 1025 K, 
then gradually drops but does not vanish even for the 
highest annealing temperature of 1400 K. The integrated 
X-ray intensity around the left position marked "twin" 
in Fig. [5fb), shown in Fig. Efb) as black squares, is also 
indicative of twinning with respect to the (111) surface 
plane. It does not only measure the surface area but also 
the average twin fraction in the film volume underneath 
the surface due to the integration depth of SXRD. The 
amount of twinning apparently does not change up to an 
annealing temperatures of 1200 K and only then grad- 
ually decreases. Roughly, both data sets agree that (i) 
twins start to decay in the temperature range of 1000 K 
to 1200 K and that (ii) even the highest annealing tem- 
perature used is not sufficient to remove them entirely. 
We attribute the shift in the onset of twin annealing to- 
wards higher temperature in the X-ray data to the higher 



pressure during Ir thin film deposition. In LEED anneal- 
ing sequences we noticed that this onset depends sensi- 
tively on the pressure during deposition. Most likely it 
is related to a higher amount of carbon containing resid- 
ual gas, which partly decomposes on the surface during 
growth. This results in incorporation of carbidic species 
into the film, which may hinder twin annealing at higher 
temperatures. 

The full dots (red) in Fig. [T^a) show the length of the 
borders separating surface areas of different stacking se- 
quence (side twin boundaries) as obtained by STM. At 
low temperature the border length is derived form deco- 
ration rows and ridges in mounds elongated along (110) 
and at high temperature from the length of the frac- 
tional steps. The border length decreases rapidly upon 
annealing at temperatures which are well below the on- 
set temperature for a decrease of the fraction of twinned 
surface area [full squares (black) Fig. Ufa)]. Apparently 
both twinned and untwinned areas become more com- 
pact (reducing side twin boundaries), and probably even 
grow in size through coalescence processes or preferential 
disappearance of small domains. Note that these pro- 
cesses take place without affecting the relative fractions 
of twinned and untwinned surface area. This is plau- 
sible, as the driving force to remove side twin bound- 
aries (separating laterally areas of different stacking) is 
much larger than the driving force to remove the coher- 
ent twin boundaries: grain boundary energies are much 
higher than stacking fault energies. The full red dots in 
Fig. E^b) measure the X-ray intensity at the central posi- 
tion marked by a vertical arrow in Fig.[6l^a), i.e. for twins 
on a {111} side-facet of a twinned crystallite, during an- 
nealing. The measured intensity displays a similar tern- 
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FIG. 7: (Color online) (a) Fraction of twinned surface area 
measured by LEED (black squares) and border length be- 
tween differently stacked surface areas measured by STM (red 
circles) for a 90 ML Ir film grown on Ir(lll) at 350 K, as a 
function of the annealing temperatures, (b) Intensities from 
the twin [black squares, measured at the left position marked 
"twin" in Fig. EJb)] and the twin of twin [red circles, mea- 
sured at the central location marked with a vertical arrow on 
Fig. [6] as a function of the annealing temperature. 



perature dependence as the border length measured by 
STM. It decreases already around 700 K, at much lower 
temperatures than the twin intensity itself. We there- 
fore conclude that the twin- of -twin intensity is linked to 
the side twin boundaries separating patches of different 
stacking. 



VI. DISCUSSION 

Below we develop a picture of the twinned film to 
account for the following facts: (i) the boundaries be- 
tween areas of different stacking sequence are parallel to 
the (110) directions; (ii) they are separated by fractional 
steps of 1/3 or 2/3 height; (hi) twins with respect to the 
{111} side facets of twinned and untwinned crystallites 
are linked to the boundaries between differently stacked 
regions; (iv) weak unexpected peaks along the CTRs par- 
allel to [111] [vertical arrows in Figs.[5fa-f)] are generated 
during the growth. 



A key to such a picture is the nature of the side 
twin boundaries. Fig. [8fa) displays an untwinned and 
a twinned crystallite bounded by (111) top facets and 
{111} as well as {100} side-facets. A cross section 
through the crystallites indicating the facets and their 
directions is shown in Fig.[8fb). Note that the directions 
are given in units of the substrate lattice vectors. There- 
fore, the [112] direction is not only normal to a {112} 
side-facet of a fee crystallite, but also to a {100} side- 
facet of a twinned crystallite. Figs. [8jc) and (d) sketch 
possible crystal plane combinations with low index planes 
forming side twin boundaries. To embed a twin crystal- 
lite into an untwinned matrix (or the opposite) one of 
the boundaries sketched in Fig. Ws) and one sketched 
in Fig. [HKd) are necessary. The principal plane combi- 
nations are of the types {112}/{112}, {001}/{122} and 
{111}/{115}. 

The {112} / {112} twin boundary is normal to the sam- 
ple surface. As known from TEM investigations of other 
fee metals, e.g. Au^ the {112} / {112} boundary leads 
to a partial dislocation of 1/2 [111] as shown in Fig.[9jb), 
i.e. the atomic positions in the twin crystal are deduced 
from the ones of the fee crystal by a rotation of 180° 
around [111], followed by a translation of 1/2 [111]. The 
addition of this translation allows a closer contact be- 
tween the two crystals along the {112} /{112} boundary. 
This would result in a fractional step height of 1 /2 mono- 
layer in contradiction to our STM observation of 1/3 and 
2/3 fractional steps. 

In a {001}/{122} boundary, every third atom of the 
{122} facet is in contact with the {001} facet and rests 
in the fourfold coordinated sites formed by the atoms of 
the {001} facet [Fig. Ms)]- This is obtained by adding 
a translation along 1/2 (110), i.e. in the surface plane, 
to the 180° rotation. The surfaces of the regular and 
twinned parts stay at the same height and thus this type 
of boundary is inconsistent with the observed partial step 
heights. 

The {111}/{115} boundary [Fig. g^d)] has one third 
of the atoms of the {115} plane in contact with the {111} 
facet. A translation of 1/3 (211), i.e. in the {111} facet, 
makes the atoms in contact sit in threefold coordinated 
fault sites of the {111} facet and increases the density 
of this boundary considerably. The same holds for the 
twin side-facets, and corresponds to a translation of 1/3 
(255), i.e. in the {115} facet. These configurations in- 
volve characteristic steps of 1/3 or 2/3 of a monolayer 
step height. 

We simulated the X-ray scattered intensity in three di- 
mensions, for a system consisting of two crystallites, one 
fee, and the other (lll)-twin, bounded by a {111}/ {115} 
twin boundary [a (111) side- facet on the fee side]. De- 
tails of the calculation are given in appendix B. Two sit- 
uations were considered: (i) the atoms of the first plane 
for each crystallite are located at positions corresponding 
to hollow sites for the (111) surface, in hep sites for the 
twinned part and in fee sites for the fee part, and the 
positions of all atoms above are derived accordingly; (ii) 




FIG. 8: (Color online) (a) Ball model of an fee and a twin 
crystallite bounded by (111) top facets and {111} as well 
as {100} side- facets, (b) Cross section depicting the facets 
and the corresponding orientations with respect to the regu- 
lar crystal matrix. On the left (right) facet of the twin (fee) 
crystallite are shown two small crystallites, one with the same 
structure as the initial twin (fee) crystallite which is displayed 
with the same colour, and one which structure is twinned with 
respect to the initial twin (fee) crystallite. The two latters 
are the twin-of-twin and side-twin and appear in light grey 
(light brown and violet). (c),(d) Models of possible side twin 
boundaries involving low index planes. 



from (i), the twinned crystallite is shifted by 1/3 [255] as 
in Fig. [9fd). The calculation along the reciprocal space 
plane denned by [111] and [112] is shown in Figs. [TOT a) 
and (b) for (i) and (ii) respectively. The most obvious 
difference between the two maps is the appearance of 
two lines, one vertical, the other horizontal (grey arrows 
in Fig. [T0|) along which constructive interferences yield 
increased intensity. The horizontal feature intersects all 
[lll]-CTRs, including also those scanned in Figs. [5fb,d) 




FIG. 9: (Color online) Ball models for high symmetry side 
twin boundaries, (a) Orientation of the crystallites: (red/left) 
twinned, (green/right) fee. (b) Model of a {112}/{112} twin 
boundary with partial dislocation 1/2 [111], (c) Model of a 
{001}/{122} boundary with a partial dislocation of 1/2 [110]. 
(d) Model of a {111}/ {115} boundary with a partial disloca- 
tion of 1/3 [255]. (e) Model of a crystal (yellow) twinned to 
a twin crystal (red) at a {111} facet, (f) Final model for the 
side twin boundary (see text). 



[marked (IV) and (V) respectively in Fig. [10]. It is worth 
noting here that the intensity is not much increased right 
at the intersection between (IV), (V), and the horizon- 
tal feature. However, if, due to some small misalignment, 
the CTRs are scanned a little off their exact position, the 
extra feature will become better pronounced. Due to lim- 
ited goniometer alignment, this is observed in Figs.[5fb,d) 
to a certain extent. Similar features were reproduced as 
well for all other cuts of the reciprocal space relevant to 
the CTRs scanned in Figs. [5fa-d). The extra features 
observed on the CTRs are thus presumably generated by 
the 1/3 [255] shifts and therefore the 1/3 and 2/3 steps. 

Fig. [9fd) provides also some clue about a possible 
origin of the strain fields accompanying the steps of frac- 
tional height apparent e.g. in Fig. [2fa). In a certain dis- 
tance from the {111}/ {115} boundary both crystallites 
must be attached to the same plane provided by the sub- 
strate - otherwise an energetically costly extended void 
would result. Attachment of both crystallites to the same 
substrate implies lattice deformations in the vicinity of 
the {111}/{115} boundary. 

We now address the side-twin and twin-of-twin peaks 
appearing along the lines parallel to [511] and [111] 
[Figs. EHa,b)]. The crystallites with densest facets corre- 
sponding to these two structures are shown in Fig. [8f b). 
The structure corresponding to the twin-of-twin is dis- 
played in Fig.[9^e). Although this accounts for the X-ray 
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FIG. 10: (Color online) Intensity map for a {111}/ {115} boundary between an untwinned and a twinned crystal of the same 
size, (a) without and (b) with a shift on the (111) side-facet resulting in a height difference of 1/3 monolayers (see text). The 
fringes are due to the finite size of the crystallites considered for the calculation. 



measurements, the (111) facet of such a doubly twinned 
domain [shown in yellow in Fig. EJe)] would have a large 
inclination of 39° [see the facet of the twin-of-twin crys- 
tallite in Fig.Ep))] degrees with respect to the (111) facet 
of twinned and fee crystallites and thus would intersect 
the thin film surface and would be seen in the STM to- 
pographs. We find no evidence for such inclined facets 
at the surface. 

A reason for this could be that the {111}/{115} 
boundary dissociates into {111}/{111} [Fig. [9je)] and 
{552}/{112} [Fig. Ep>)] boundaries, as sketched in 
Fig. MS)- This is well supported by the literature, 
where the instability and dissociation of the {111}/ {115} 
boundary is reported 2 ^. Note that such a dissociation 
still assumes that the twin and fee domains are trans- 
lated with respect to each other with a 1/3 [255] shift 
which is responsible for the extra feature discussed ear- 
lier. Energetically, the dissociation proposed is consistent 
with the available literature data for other fee systems. 
For instance Wolf et al 28 calculated for Cu the energy 
of a {111}/{115} boundary to be 0.61 Jm -2 , i.e. larger 
than the sum of the energies for a {111}/{111} bound- 
ary (0.01 Jm" 2 ) and for a {112}/{552} boundary (0.26 
Jm" 2 ). 

The bottom facet of the light grey (yellow) crystallite 
visible in Fig. [9jf) is a {113}. If laterally extended, it 
could be matched by a {335} facet growing out of the 
(111) substrate to form a high symmetry boundary. The 
{113}/{335} boundary was estimated to have an energy 
about four to 20 times larger than that of a {111}/{111} 
boundary- 2 ^—. In other words the {111}/ {115} bound- 
ary dissociation also implies the formation of a costly 
{113}/{335} replacing a {111}/{111} boundary parallel 
to the substrate surface. Yet, the energetic cost is lim- 
ited due to the short length of the {113}/{335} boundary 
and appears compensated by the formation of the other 
{111}/{111} boundary (the side one). 

Although the dissociated {111}/ {115} structure is en- 
ergetically more favorable than the undissociated one, 
it is considerably less favorable than the coherent (111) 
twin boundary. Therefore upon annealing the driving 



force to remove the dissociated {111}/ {115} boundaries 
is much larger than the one for the removal of the coher- 
ent twin planes. Thus first the length and extension of 
the {111}/ {115} boundaries diminishes through coars- 
ening as observed by STM and SXRD, and only then, 
at higher temperature the twins are eventually removed. 
We speculate that the removal of the twins taking place 
at about half of the melting temperature in itself is linked 
to the onset of bulk diffusion through the onset of vacancy 
generation. 

Considering one twin crystallite, the discussion was fo- 
cussed on two of its boundaries with the remaining crys- 
tal: the twin boundary in the (111) plane, and one of the 
side boundaries. However, minimization of the system's 
energy may require maximum contact of the twin grain to 
its three dimensional matrix, noticeably through contact 
planes such as the one for the dissociated {111}/{115} 
boundary. This may for instance imply strain fields in 
the grain and surrounding. Such an issue is beyond the 
scope of the article. 

VII. CONCLUSION & SUMMARY 

Homoepitaxial thin Ir films growth at room temper- 
ature on Ir(lll) is characterized by a proliferation of 
twins starting from 10 ML, leading to an areal frac- 
tion of twin approaching 70 % for about 70 ML, as 
determined by SXRD and STM. Quantitative analysis 
of SXRD data from relatively thick films (several 10 
ML) was accordingly developped and the twin/fee com- 
position of the film was obtained. The boundaries be- 
tween different stacking areas are identified to consist 
of {111}/ {115} boundaries probably dissociated into co- 
herent {111}/{111} and {552}/{112} boundaries. This 
transformation takes place on the {111} side-facets. It 
involves the formation of a volume twinned with respect 
to twin crystallites. The proposed structure of the dis- 
sociated boundaries is consistent with the observed ad- 
ditional peaks in the the SXRD data and the STM ob- 
servation of fractional step heights of 1/3 and 2/3 of a 
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monolayer along (110). 

During annealing of the Ir thin films first the area 
of the boundaries separating different stacking dimin- 
ishes through coarsening at temperatures in the range of 
800-1000 K. Only at much higher temperatures beyond 
1200 K also the twins themselves heal. The two step an- 
nealing process is traced back to a difference in the driv- 
ing force for healing: the boundaries between different 
stacking domains are energetically more costly than the 
coherent twin boundaries. 
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Appendix A: Details of the simulation of the X-ray 
intensity scattered along CTRs 
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34% 




equivalent 
stacking 




The X-ray attenuation length was taken into account in 
the calculation through an imaginary part3 for the scat- 
tering vector resulting form the absorption index which 
adds to the refraction index of X-rays. 

The calculation includes the contributions from the 
bulk and the thin film on top, as well as the interfer- 
ences between them. Since the lateral extension of the 
twin and fee regions, as assessed by STM, is much smaller 
than the X-ray beam coherent length, their interferences 
are included as well. For the calculation, 10 vertical do- 
mains (columns) are considered. The twin nucleation at 
different stages of the thin film growth is accounted for 
as follows. In each of the 10 domains one stacking fault, 
initiating the formation of a twin, is formed at a given 
height. If i is the index of the domain, i G [1, 10], and t is 
the film thickness, the stacking fault is formed at a height 
i x ip (^) , ip being the integer part function. Below this 
height the domain consists in a pure fee crystal, above it 
consists in a fraction / of twin and the remaining of fee 
(see Fig. [11] for an illustration with three domains). If all 
domains were given the same weight in the summation 
of their structure factors, the overall fee/ twin fraction 
would increase linearly along the growth direction (in- 
clined dashed line in Fig. [11]). So that this increase can 
be slower or faster, the domains may all have a different 
weight, i.e. a different lateral size. Instead of using inde- 
pendent weights, which would imply 10 parameters to be 
refined along the fit procedure, we impose that the weight 
Wi of a domain depends on its index, i.e. Wi = i 1_r? , rj 
being the only parameter to be refined. This is exempli- 
fied in Fig. HI1 

From this model one may deduce the equivalent 
twin/fee fraction at a given height in the thin film by 
adding the fee and twin contributions from each of the 
10 domains at this height. For the atomic layer la- 
belled n in the film (n = 1 for the layer in contact 



FIG. 11: (Color online) Schematics of the model used for the 
X-ray scattering simulations. The twinned film is sliced in 10 
columns (for the sake of clarity this is exemplified here with 
three) . Each of the columns comprises fee and twinned stack- 
ing, in light gray (green) and dark gray (red) respectively. 
As the index of the column increases, the first twinned re- 
gions (coexisting with fee ones in a f/1 — f proportion with 
/ = yjjjj) are found higher in the column. Each column is 
given a different lateral size, which is accounted for by the ex- 
ponent rj = 0.7 (see text). The equivalent stacking along the 
growth direction is built by arranging twin regions from each 
column side by side, and the same for fee regions. The dis- 
tinct lateral sizes of each column result result in a non- linear 
(dashed line) evolution of the twin fraction as a function of 
the height in the thin film. 



with the substrate), the fec/twin fraction is constant 
if n e [ip(±)x(j-l),ip(±)xj] with j e [1,10] 
since new twinned domains nucleate only at multiples of 
ip (^) . So for n e [ip (^) x (j - 1) , ip (^) x j] , the 
twin fraction is 



and the fee fraction is 



F n = 



(i-/)n=i^-"+£r= J+ ifc i - r ' 



(Al) 



(A2) 



with the normalization factor 



s = (i-jOELi^ 1 -" 

= nil fc 1 -" 



I-77 



(A3) 

Figure 03(h) was derived from these equations. 
Note that the choice of 10 vertical domains (and thus of 
stacking fault nucleation at multiples of ip(t/ 10) in each 
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domain) is a compromise between computation time and 
fine description of the system. The use of more than 10 
domains would have been prohibitively time consumming 
along with the fitting procedure, which implies calculat- 
ing the intensity a large number of times. 

The roughness is included following a /3-model, namely 
the film does not end abruptly but rather a number of 
layers are added on top which occupancy decays with the 
index number n of the layer following /3 n , f3 G ]0, 1[ be- 
ing the roughness parameter. The j3 parameter is chosen 
the same for the twin and fee regions. f3 may be con- 
verting into an rms roughness by calculating jz^di r , dj r 
being the lattice spacing perpendicular to the surface. It 
should be noted that we found that translations of some 
regions in the crystal ensure the crystallite contacts at 
twin boundaries. These translations have a component 
perpendicular to the substrate corresponding to shifts of 
a fraction of the lattice spacing. In addition, the deco- 
ration rowsii also are regions of the crystal into which 
the atomic positions are displaced perpendicular to the 
surface with respect to the atomic positions in fee and 
twin regions. Such contributions (fractionnal shifts and 
decoration rows) of a non integer multiple of dj r are in- 
cluded in the rms roughness evaluation by STM, and not 
in the model for the simulation of the CTRs scattered 
intensity. This makes the direct comparison between the 
values derived by the two techniques very unreliable. 

Four parameters are refined during a fit: the twin frac- 
tion /, the exponent r] for the domains weights, the f3 



roughness parameter, plus a scale factor between the data 
and the simulation. The fit is performed after the data 
is corrected taking into account the Z-axis geometry^ of 
the goniometer which is relevant here. 



Appendix B: Details of the simulation of the 
contribution of {111}/ {115} twin boundaries to the 
X-ray scattered intensity 

For simplicity, the calculation was here performed 
without taking into account any attenuation length for 
the X-rays. The structure factor was computed for a 
three dimensional crystal including a fee part and the 
remaining twinned along [111]. The side of the fee part 
closeby the twin one ends with a {115} facet and the cor- 
responding facet for the twin is a {111}. Two situation 
were considered, with or without a 1/3 [255] translation 
thanks to which the fee and twin parts may be brought 
in contact. For the two-dimensional maps of the scat- 
tered intensity presented in this article, the dimensions 
of the overall crystal (twin+fee) were 130 A along [111], 
and 110 A along both [211] and [112]. The generality of 
the qualitative conclusions drawn in the text do not de- 
pend on the relative fraction of twin and fee or size of 
the crystal. 

Since the crystal's ends are perfectly abrupt, finite size 
fringes develop in the intensity maps. 
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